BOUNDARIES OF DISK-LIKE SELF-AFFINE TILES 



KING-SHUN LEUNG AND JUN JASON LUO 

CN . Abstract. Let T := T(A, V) be a disk-like self-affine tile generated by an integral 

expanding matrix A and a consecutive collinear digit set £>, and let f(x) = x 2 + 
^vq ' px + q be the characteristic polynomial of A. By investigating the neighbors of 

(— i ■ T, we construct a neighbor graph and show that any point on the boundary dT 

^ \ can be identified with a one-sided infinite path of the graph. Also we derive 

certain equivalent conditions for the pair (A, T>) to be a number system. From the 
neighbor graph we can deduce a family of graph-directed sets, which gives a second 
description of dT. When A is a similarity, we obtain dim jj(dT) — 2 log pj log \q\ 
■ where p is the largest positive zero of the cubic polynomial x 3 — (\p\ — l)x 2 — (\q\ — 

\p\)x- \q\- 

-I— > 

1. Introduction 

Let M n (Z) denote the set ofnxn matrices with entries in Z and let A e M n (Z) 



X 



be expanding (i.e., all eigenvalues of A have moduli > 1). Assume |det(A)| = \q\ 



00 , 

and P = {0, di, . . . , C Z n with |g| distinct vectors. We call V a c?^«t set and 

(A, P) a self-affine pair. It is well-known that there exists a unique self-affine set 
O ! T:=T(A,V) [LWT] satisfying 



T = A- 1 (T + P) = |gA-^:rf Ji Gl? 

If T has non-void interior, then there exists a subset J7" C Z n such that 

T + J = M n and (T + t)° n (T + t')° = 0, t ^ t', t, t' e J", 

thus T is called a self-affine tile and J7" a tiling set. T + J is called a tiling of M n , 
and a lattice tiling if J is a lattice |LW3j . 

The topological properties of self-affine tiles, such as connectedness or disk-likeness 
(i.e., homeomorphic to the closed unit disk), have attracted a lot of interest. A 
systematical study on connectedness was due to Kirat and Lau |KL] . they mainly 
concerned a class of self-affine tiles T generated by the consecutive collinear (CC) 
digit sets V := T>(v, \q\) = {0, 1, . . . , \q\ — l}v, v G Z n \ {0} via the algebraic property 
of the characteristic polynomial of the matrix A. Other considerations on non- 
consecutive collinear or non-linear digit sets can be found in |DLj . |LLu] . The 
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question of disk-likeness was first investigated by Bandt and Gelbrich |BG] for self- 
affine tiles in M 2 with | det(A)| = 2 or 3. They also observed that the characteristic 
polynomial of an expanding matrix A e Mii^) is of the form: 

f(x) = x 2 + px + q, with \p\ < q, if q > 2; \p\ < \q + 2|, if q < —2. 

By studying the neighborhood structure of T, Bandt and Wang |BW] proved that a 
tile T with no more than six neighbors is disk-like if and only if T is connected. A 
translation of the tile T+£, £ e J is called a neighbor oi T if Tn(T+£) ^ 0. Making 
use of this criterion, Leung and Lau |LL] then gave a complete characterization of 
the disk-likeness of self-affine tiles with CC digit sets. 

Theorem 1.1. f |LL] j Let A e M 2 (Z) be an expanding matrix with characteristic 
polynomial f(x) = x 2 + px + q. Then for any CC digit set TJ(v, \q\) in 7? such that 
v, Av are linearly independent, T is a disk-like tile if and only if 2\p\ < \q + 2|. 
Moreover, when p = 0, T is a square tile; when p ^ 0, T is a hexagonal tile. 

Recently, Akiyama and Loridant ( [AL] . |AL2j ) proposed a method to parameterize 
the boundary of a self-affine tile, and reproved Theorem 11.11 by showing that the 
boundary of T is a simple closed curve. In the present paper, we go further to explore 
the structure of the boundary of the T defined in Theorem ll.il For convenience, we 
call such T a CC tile. If it is also disk-like, we call it a disk-like CC tile. 

As the first main result of this paper, we describe the boundary dT of a disk-like 
CC tile T in two equivalent ways. First we define a labeled directed graph (called 
the neighbor graph of T) with V = {£ : T + i is a neighbor of T} as the set of 
vertices and AT> the set of edge labels, such that dT is identified as the union of all 
one-sided infinite paths of this graph. Hence dT determines a sofic system [Fij. The 
neighbor graph is associated with a contact matrix M whose entries are the edge 
labels. Then we provide some equivalent conditions for the self-affine pair (A, T>) to 
be a number system |MTT] . On the other hand, dT can also be expressed as a family 
of graph-directed sets |SWj by using an algorithm (Proposition 13. II) . We prove that 
the associated graph-directed IFS satisfies the open set condition. From this we can 
deduce our second main result, an expression for the Hausdorff dimension of dT. 
The derivation involves finding the Perron-Frobenius eigenvalue of M. When A is 
a similarity, we obtain the following. 

Theorem 1.2. Let A e M 2 (Z) be an expanding similarity with characteristic poly- 
nomial f(x) = x 2 + px + q and T = T(A, T>) be a disk-like CC tile. Then 

dimH (ST) = 

log \q\ 

where p is the largest positive zero of the cubic polynomial x 3 — (\p\ — l)x 2 — (\q\ — 
\p\)x- \q\. 

The rest of the paper is organized as follows: in Section 2, we construct a neighbor 
graph and identify dT with a sofic system determined by this graph, moreover we 
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derive some equivalent conditions for the self-affine pair (A,T>) to be a number 
system; in Section 3, we consider dT as graph-directed sets and prove Theorem 11.21 
by checking the spectral radius of the contact matrix. Finally all neighbor graphs, 
graph-directed sets and contact matrices corresponding to different characteristic 
polynomials f(x) are listed in Appendices A, B and C for easy reference. 

2. Boundary of T as a Sofic System 

We first introduce some terminology of symbolic dynamics from [LMj . Let Q = 
Q(V,£) be a directed graph where V is the set of vertices and £ the set of edges. 
Let A be a finite set (called alphabet). If there exists a mapping (called labeling) 
C : £ — >■ A, then the ordered pair G = (G,C) is called a labeled directed graph. All 
the infinite paths £ = eie 2 e 3 ... on ^ constitute the so-called edge shift "Kg. Define 
the label of the path £ by 

C oo {0:=£(e 1 )C{e 2 )jr{e 3 )...eA n . 

Here *4. N is called the full shift of A. The set of all such labels is denoted by 

X G = {xeA N :x = £«,(£) for some £ G X g } . 

Any subset of ^4 N which can be defined by a labeled directed graph as above, is 
called a sofic shift or sofic system ([Fij, [LMj ). Weiss |We] coined the term sofic 
which is derived from the Hebrew word for finite |LMj . 

Let D = {0, 1, . . . , |g| - 1} and the difference set AD := D-D, then the CC digit 
set T> = Dv and AT> := D — T> = ADv. Without loss of generality, we can assume 
the digit set V is primitive [LVV2] (i.e., Z 2 = {jv + 5Av : 7, 5 G Z}) and T + Z 2 is a 
lattice tiling. It is easy to see that T + i is a neighbor of T if and only if 

00 

I = 1V + SAv = biA~*v G T - T, where b { G AD. 
i=i 

The following is a neighbor-generating formula which plays a key role in constructing 
the labeled directed graph for the boundary. 

Lemma 2.1. ( \LL\ ) Suppose T + i is a neighbor of T with I = jv + SAv = 
hA~ l v, then we get another neighbor T+£' satisfying £' = At—b x v = j'v+5'Av 
with 7' = — (qS + bi) and 5' = 7 — pS. 

Inductively, we can construct a sequence of neighbors: {T + i n }^Lo where £0 = £ 
and £ n+ i = Al n - b n+1 v. 

Let T be a disk-like CC tile and V = {£ : T + £ is a neighbor of T}. Let Tg = 
T H (T + £). Then the boundary of T can be written as 

(2.1) dT=\jT e . 

lev 

Define an edge set £ := {e = (£, £') : £, £' G V and = v4£ — &ii> for some bi G 

AD} and a labeling £ : £ — > A by £(e) = b\ where A = AD. Then by the definition 
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above, G = (Q, C) is a labeled directed graph and it determines a sofic shift. We 
call G the neighbor graph of T. 

Proposition 2.2. Let G be the neighbor graph of a CC disk-like tile T. If x = 
YuiLi a iA~ l v = £ + YuiLi a \A~ % v G Ti where aj, a\ G D, then {bi := Oj — a' i \°l l is the 
sequence of labeling of the edges of an infinite path starting at £ (or simply called a 
label sequence starting at I). Conversely, any label sequence {h}^ (with bi G AD) 
starting at I defines a set 

OO oo 

{x : x = OiA~ % v = i + a[A~ l v, — a[ = bi, a i; o! i G D for % = 1, 2, . . . } 

i=l i=l 

of boundary points of T . 

Proof. Since I = J^i=i biA~ l v with 6j = — a[, by Lemma [2.1[ we have a sequence 
of neighbors {T + £ n }^o where £ = £ an d C+i = A£ n ~ b n +iv, hence {b i } c *L 1 is a 
label sequence starting at £ by the definition. 

Conversely, if £ = Yli=i biA~~ l v where bi G AD, then bi = — for a^, a[ E D 
and £ = 5^=i( a » — a^)A~V It follows that 

oo oo 

(2.2) x = a i A ~ iv = i + Yl a 'i A ~ iv ^ T H (T + £) = Ti. 

i=l i=l 

□ 

We can verify whether the origin is a boundary point of T by the following 
method. 

Corollary 2.3. G dT if and only if there exists an infinite path in G with all edge 
labels either non-positive or non-negative. 

Proof. Suppose G T fl (T + €) for some neighbor T + £. Putting cij = for all i 
into (12. 2p . we have 

oo 

i = £(-aJ)A-s/. 

i=l 

Since a\ G D, the label sequence {fej = —a' i }°^ 1 starting at £ has all labels non- 
positive. Similarly {b'i = ai}^ is a sequence starting at — £ with all labels non- 
negative. By reversing the argument, we can prove the converse. □ 

In fact, we can determine the neighbor graph G for any disk-like CC tile T. Let 
us take the case of f{x) = x 2 +px + q, p, q > 2 (excluding p = q = 2) as an example. 
By Theorem 11.11 T is a hexagonal tile with six neighbors |LL] and 

(2.3) V = {±v, ±(Av + {p-l)v), ±(Av+pv)}. 

In view of the definition of S , if i — v we take b\ = —p and £' = Av + pv or 
bi = —(p — 1) and £' = Av + (p — l)v; if £ = Av + pv, using f(A)v = 0, we have 
bi = — (g— 1) and £' = —v. Proceeding similarly with all i, we obtain Table [U Then 
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£ 


h 


f 


V 


-(p-1) 


+ (p — l)v 




-p 


Av + pt> 


Av + (p — l)v 


-(q-p) 


—Av — pv 




-(q-p+1) 


—Av — (p — l)v 


Av + pv 




—v 


—v 


p — 1 


—Av — (p + l)v 




P 


—Av — pv 


—Av — (p — l)v 


g-p 


Av + pv 




g - p + 1 


Av + (p — l)f 


—Av — pv 


g-1 


f 



Table 1. Relation among all neighbors of T associated with f(x) = 
x 2 + px + q, p, q > 2, 2p < q + 2 (excluding p = q = 2). 







Figure 1. The neighbor graph of T associated with f(x) = x 2 +px + 
q, p, q > 2, 2p < q + 2, (excluding q = p = 2). 

we establish the neighbor graph (Figured]). The neighbor graphs corresponding to 
other f(x) are given in Appendix A. 

Following jLWlj . we let T> Ajk = {^to a i Aiy '■ a i e D }i AT> A,k = ^A,k - ^A,k = 
{Eto M*« : 6< e AD} and P Aj00 = U^P^. 

Proposition 2.4. Lei T be a disk-like CC tile and T + £ a neighbor. Then £ = 
Y^i = Q°iA % v E AT>A,k+i for some k E Z ™tt G { — 1, 1}, h E AD for < i < k. 
When f(x) = x 2 ± 2x + 2, k = 3; and k = 1 otherwise. 
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Proof. It follows from ( I2.3P that £ G AT>a,2 excluding the case of f(x) = x 2 ±22 + 2. 
For f(x) = x 2 ± 2x + 2, we have Av ± 2v = A 3 v ± A 2 v + Av e AVa,a by using 
(A T I)f(A) = 0. □ 

A more desirable property is that any leZ 2 can be expressed as £ = Yli=o o,iA % v G 
fc + i (instead of AZ^^+i) for some fceZ with = 1 and aj G -D where < z < k. 
But this is not always the case. This property is closely related to a number system 
defined below (see also |MTTj ). 

Definition 2.5. Let A G M^i^) be expanding and V be a CC digit set. The self- 
affine pair (A,D) is said to be a number system if for any £ G Z 2 ; it has a unique 
representation £ = X]i=o^X' w ^ v '% ^ ^ 

For convenience, we sometimes write a point of the form x = ^A^v G M 2 
as radix expansion: 0.010203 . . • An overbar denotes repeating digits as in 0.12301 = 
0.12301301301 . . . Likewise, o_20_iao- 010203 . . . represents a point a_2A 2 v + a-iAv + 
a o v + 2~2i^=i OjiA~ l v. Note that shifting a radix place to the left means multiplying A 
to x. When x is on the boundary of T, the radix expansion of x is not unique. Now 
we give some equivalent conditions for the self-affine pair (A, T>) to be a number 
system. 

Theorem 2.6. Let T = T(A,V) be a disk-like CC tile. Then the following are 
equivalent: 

(i) (A, V) is a number system. 
(11) G T°. 

(Hi) f(x) = x 2 + px + q with —l<p and q > 2. 

(iv) For all neighbors T + £, £ = Yli=o a iA t v G £>A,fc+i for some k G Z with = 1 
and a,i G D where < i < k. 

Proof. (i)=Kii) Suppose £ T°. Then G In (T + £) for some £ G Z 2 \ {0}. Since 
(A, X>) is a number system, £ = Y^.=-k a iA~ l v with Oj G -D and a_k > 0. Hence 
= a_fca_(fc_i) . . . a_iOo.oi0203 .... Shifting the radix point k places to the left, we 
get = a_fc.a_(fc_i) . . . a_iaodia2Ci3 .... That means T + a_fct> is a neighbor of T. By 
Proposition I2.4[ a -k — 1- Hence corresponds to an infinite path starting at v with 
non-positive labels 6j = — o^. But by checking all the neighbor graphs in Appendix 
A, we find no such path. 

(ii) =>(i) It suffices to show that Z 2 C T>a,oo- By the lattice tiling property, 
is the only lattice point in T, i.e., Z 2 n T = {0}. It follows that Z 2 n A n T = 
l^JiZo A % V = T>A, n fo r n > 1. If i G Z 2 , there exists a large integer n such that 
£ G A n T as G T°, then £ G P A ,n C V Aoo . 

(ii)-^(iii) By inspecting all neighbor graphs in Appendix A, we find that in each 
graph corresponding to f(x) = x 2 + px + q with — 1 < p and q > 2, there exists 
no infinite path with edge labels either all non-positive or all non-negative, hence 
G T° by Corollary 12.31 In every other case, there always exists such a path. All 
these paths are listed in Table [2j 
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Six) 


Neighbor 


Path 


x 2 - 2x + 2 


— v 
—Av + u 


(-1) 
T 


x 2 — px — q 


V 

—v 


p{q - 1) 
(-p)l-(q-l)} 


x 2 + px — q 


Av + (p+ l)v 
—Av — (p + l)t; 


(q-p-l) 
\-(q-p- 1)1 

L \l Mr /J 


x 2 — q 


f 

— V 
Av + u 
— Av — v 


O(g-l) 
Of— fo — 1)1 

(9-1) 
[-(9-1)] 


x 2 — px + q 


Av — (p — l)v 
—Av + (jp — l)v 


[-(g-p + i)] 
(g-P + i) 



Table 2. Infinite paths representing a boundary point 0. 



(iii) =>(iv) Let /(x) be one of the cases: x 2 + q, x 2 + x + q, x 2 + px + q (p > 
2, excluding p = q = 2), x 2 + 2x + 2, x 2 — x + g, where p > and q > 2. In each 
case, we can rewrite their neighbors as the desired form in (iv). By using = f(A)v, 
= (A - I)f(A)v, = (A + I)f(A)v, we have 

Case (1) f(x) = x 2 + q. Av — v = A 2 v + Av + (q — l)v, —v = A 2 v + (q — 
l)v, -Av = A 3 v + (q - l)Av, -Av + v = A 3 v + (q - l)Av + v, -Av - v = 
A 3 v + A 2 v + (q- l)Av + (q - l)v. 

Case (2) f(x) = x 2 + x + q. -v = A 2 v + Av + (q - l)v, -Av = A 3 v + A 2 v + 
(q — l)Av, —Av — v = A 2 v + (q — l)v. 

Case (3) f(x) = x 2 +px+q (p > 2). —v = A 2 v+pAv+(q— l)v , —Av — (p—l)v = 
A 2 v + (p — l)Av + (q — p + l)v , —Av — pv = A 2 v + (p — l)Av + (q — p)v. 

Case (4) f(x) = x 2 + 2x + 2. Av + 2v = A 3 v + A 2 v + Av, -v = A 4 v + A 3 v + 
A 2 v + v, -Av - v = A 2 v + Av + v, -Av - 2v = A 2 v + Av. 

Case (5) f(x) = x 2 — x + q. Av — v = A 2 v + (q — l)v, —v = A 3 v + (q— l)Av + (q — 
l)v, -Av = A i v + (q-l)A 2 v + (q-l)Av , -Av+v = A 4 v + (q-l)A 2 v + (q-l)Av+v. 

(iv) ^>(ii) Suppose ^ T°. By the same argument as in the proof of (i)=^(ii) 
above, there should be an infinite path in the neighbor graph starting at v with edge 
labels all non-positive. But we find no such paths by inspecting all the neighbor 
graphs in Appendix A. □ 

Remark 2.7. Gilbert [GiJ obtained some related results in the context of quadratic 
number fields. We conjecture that Theorem \2.6i can be extended to non-disk-like tiles. 
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3. Boundary of T as Graph-directed Sets 

For a directed graph Q = Q(V,S) where V = {vi,v 2 , ■ ■ ■ ,v m }, we write £ it j for 
the set of edges from vertex Vi to vertex Vj, and we add a contraction mapping 
F e : IR 2 — > IR 2 for each edge e G £. Then the family of contractions {F e : e G £} is 
called a graph- directed iterated function system (graph- directed IFS) and it is well- 
known that there exists a unique family of non-empty compact subsets E\, . . . , E m 
of R 2 ([FES, pW] ) such that 

m 

(3.1) Efc = (J U F ^ E i)- 

j=l eeSij 

We call {E\, . . . , £ m } a family of graph- directed sets. The graph-directed IFS {F e : 
e G £} is said to satisfy the open set condition ( OSC) if there exist a family of open 
sets {0\, . . . , O m } such that 

m 

(3.2) Oi d |J |J F e {Oj) fori = l,2,...,m 

i=i ee£ ii3 

with disjoint unions, i.e., F e (Oj) fl F e i(Oji) = whenever (e, j) 7^ (e',f). 

In this section, we will show that the relationship of {T^gy where V = {£ : 
T + £ is a neighbor of T} as in Section 2 can be expressed as a family of graph- 
directed sets by means of the following proposition, which applies to non-disk-like 
T as well. 

Proposition 3.1. Let £ = + 5Av, £' = j'v + 5'Av G V such that £' = Al - b x v 
for some b\ G AD, then 

A-lfrr I • ^ r- T t 11 J {&1, &1 + 1, • • • , ? ~ 1} */ &1>0; 

(T?+F)CT< /orflM := i {0,l 5 ...,g-l + &1 } if bl <0. 

Moreover, 

e'eBijei bl 

where B e := {£" G V : £' ' = Al - b[v for some b[ G AD}. 

Proof. When b% > 0, if x G Ty then the radix expansion is 

x = 0.cic 2 c 3 . . . = S'j' ' .d x d 2 d z .... 

It follows from Lemma [2.11 and = A~ 1 f(A)v that 

A~ x x + (61 + k)A~ l v = 0.(6i + k)cic 2 c z ... = 8j.kc[c 2 c' 3 ... eT e 

for k = 0,1, . . . ,q — 1 — b\. The case when b\ < can be proved similarly. 
For the second part, we only need to show 



Let y = O.a^a^ . . . = Sj.a^ai^a^ . . . G T e . It follows that Ay — a x v = 0.a 2 a 3 a 4 . . . = 
#7(0^ — ai).a' 2 a' 3 . . . G Tp, where £' = A£—(ai — a' 1 )v. This implies y G A^ l {T^ +a\v). 
By definition, we see that a\ G h 1 for b% — a\ — a' v □ 



Therefore, dT is the union of a family of graph-directed sets. We give the case 
f(x) = x 2 + px + q (p, q > 2 excluding p = q = 2) as an example in the following. 
All the other cases are given in Appendix B. 

Example 3.2. Consider the case f(x) = x 2 +px + q (p, q > 2, excluding p = q = 2). 
When £ = v, from Table[J\ we have = B v = {Av + pv,Av + (p — 1)^}. When 
£' = Av + pv , bi = —p and I_ p = {0, 1, 2, . . . , q — 1 — p}; when £' = Av + (p — l)v, 
bi = —{p — 1) and I-(p-i) = {0, 1, 2, . . . , q — p}. Thus by Proposition ^. 1\ the first 
set equation comes out. Similarly the other five can be deduced. For simplicity 
we let ui = v , u 2 = Av + (p — l)v, u 3 = Av + pv. Then the graph- directed sets 
T± Ul , T± U2 , T± U3 , representing dT satisfy 



q-p q-p-1 

+ jv) U [J (T U3 + jv) 

j=0 j=0 
p-2 p-1 

"3 3 V ) 

j=0 j=0 

T 

q-l q-l 

Us + j V ) 

j=p-l j=p 

q-l q-l 

|J (T U2 +jv)U |J {T U3 +jv) 

j=q-p+l j=q-p 

T U1 + l ) v 



Theorem 3.3. The graph- directed IFS representing the boundary of a disk-like CC 
tile T satisfies the OSC. 

Proof. Replacing Ti by (T + £)°, we can check that the graph-directed IFS repre- 
senting dT satisfies the OSC. We illustrate the idea by proving the case f(x) = 

x 2 + px + q (p > 2, q > 2, excluding p = q = 2). In view of Example 13.21 we need 
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AT Ul 



AT 



AT 



"3 



AT- Ul 



AT- 



"2 



AT_ 



"3 



to show 

q-p q-p-l 

A{T + Ul )° D \J((T + u 2 y+jv)U |J ((T + u 3 ) +jv) 

j=Q j=0 
p-2 p-1 

A(T + u 2 )° D \J((T-u 2 y+3v)u\J((T-u 3 )°+3v) 

j=0 j=0 

a{t+u 3 )° d (r- Ul y 

q-l q-l 

A(T-u x y D |J ((T-u 2 y+jv)u\J((T-u 3 y+jv) 
j=p-l j=p 

q-l q-l 

a(t-u 2 )° d u ((t+ M2 )°+^)u u ((T+^r+jV) 

j=?-p+i j=q-p 
A(T-u 3 y D (T + Ul )° + ( g -l> 

with disjoint unions. Since T is a CC tile, it follows that 

q-l q-l 

(3.3) AT° d |J (T + jv)° = |J (T° + jv) 

i=o i=o 

with disjoint union. By using (13. 3p and = = A 2 v + pAv + qv extensively, 

we prove the first two set inequalities in the following. The remaining four can be 
verified similarly. 

For j = 0,1, ...,q-p, 

(T + u 2 )° + jv = T° + (p - 1 + j)v + AvC A(T + Ul )°. 
For j = 0, 1, . . .,q-p - 1, 

(T + u 3 )° + jv = T° + (p + j)v + AvC A(T + Ul )° . 
For j = 0,1,..., p-2, 

(T-u 2 )°+jv = T° + (j -p + l)v- Av 

= T° + (q + j -p+ l)v + A 2 v + (p - l)Av 
C A(T + u 2 )°. 

For j = 0, 1, . . .,p - 1, 

(T — w 3 )° + jv = T° + (j — p)u — Aw 

= T° + (q + j - p)v + A 2 v + (p - l)Av 
C A(T + m 2 )°. 

By the same way, all the other cases follow and hence the theorem is proved. □ 

The Hausdorff dimension (dim#) (see e.g., [Faj . |Fa2j ) is the most common and 
important dimension in fractal geometry. Since a disk-like tile has positive Lebesgue 
measure, its Hausdorff dimension is simply the dimension of the underlying space 
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R 2 . A more interesting problem is to find the Hausdorff dimension of the boundary 
of the T. In the rest of this section, we find the exact value of dinif/(<9T) for certain 
particular cases that A is a similarity. We state the simplest one first. 



Proposition 3.4. Let A e M 2 (Z) be expanding with characteristic polynomial 
f\x) =x 2 + q (\q\ > 2) and T(A,V) a disk-like CC tile. Then dim^(<9T) = 1. 

Proof. By Theorem II .![ T is a square tile (parallelogram). Hence dim#(<9T) = 1. □ 

Geometrically, a similarity is a multiple of either a reflection or a rotation. We 
call the former a scaled reflection and the latter a scaled rotation] algebraically, a 
similarity is a multiple of an orthogonal matrix. The case that A is a scaled reflection 
is solved already as its characteristic polynomial is of the form f(x) = x 2 — q (q > 0). 
So we focus our attention on those A that are scaled rotations. 

Lemma 3.5. Let A be a scaled rotation. Then its characteristic polynomial has 
positive constant term and A has either two distinct non-real eigenvalues or two 
equal real eigenvalues. 

Proof. Let A = ( rcos ^ rsinO \ characteristic polynomial is given by x 2 — 



rsin^ rcos# 

2rcos6 x + r 2 . It has two equal real zeros when 9 = or it and two distinct non-real 
zeros otherwise. □ 

We shall make use of the following well-known dimension formula (see also |DKV] , 
[Fa2] . [HLR] . [EN] . [MWj . [SW] ) to compute the Hausdorff dimension of dT . 

Theorem 3.6. If A is a similarity with |det(v4)| = \q\ > 2, then the Hausdorff 
dimension of dT is given by 

(3.4) dim^(<9T) = log p{M)/ log r = 2 log p{M)/ log |g|, 

where p(M) denotes the spectral radius of the contact matrix M and r = |g| 1,/2 is 
the expansion ratio of A. 



Remark 3.7. In |MW] , a separation condition is required for the validity of flff.^p . 
It is relaxed to the OSC in ^2]. It is shown in jDKVj . jHTR] . [LN] and [SW] tfiat 
jjj3.4\) is also valid with a weaker assumption of the OSC. 



We first find the contact matrix (or substitution matrix) M for T. Let i,£',bi 

and Bi be defined as in Proposition 13.11 The entry Mm is defined by Mu> : = 

#((A£ + V) (1 (£' + V)) . Since V is a CC digit set, it can be shown easily that 

Mu' = = q— \bx\ where 1^ is as in Proposition 13.11 Recall that b\ is the label 

of the edge starting at I and ending at £'. Hence we obtain the contact matrix M 

of T from its neighbor graph with different edge labels (i.e., replace b% by q — |&i|). 

There is a one-to-one correspondence between the contact matrix and the neighbor 

graph. For example, the contact matrix for the case f(x) = x 2 + px + q (p,q > 

2, 2p < q + 2 excluding p = q = 2) can be found in Table [3J and the related 
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V 


/J 1 1 1 1 ry\ \ 111 

-tiv -\- {p ±)u 


i±u ^ pu 


— V 


/\ i ! I n I 1 o i 

— VT. L' — — _L J U 


/l 1 1 I'll i 
— -r± U — pv 


V 





q-p + 1 


q-p 











Av + (p — i)v 














p — 1 


P 


Av + pv 











1 








—v 














q-p+1 


q-p 


—Av — (p — l)v 





p — 1 


P 











—Av — pv 


1 


















Table 3. The contact matrix of T associated with f(x) = x 2 + px + 
g, p, q > 2, 2p < q + 2 (excluding p = q = 2). 



neighbor graph is shown by Figure [TJ The contact matrices for the other cases are 
given in Appendix C. 

If M is irreducible (i.e., for each entry My, there exists an integer n > such that 
(M n )ij > 0), then the spectral radius p(M) = Xm where Am is the Perron- Frobenius 
eigenvalue of M as stated in the following simplified version of the Perron- Frobenius 
Theorem (see |Ga] or [Se] for details). 

Theorem 3.8. Let M be an irreducible non-negative matrix. Then there exists a 
positive eigenvalue Xm such that Xm > M for all eigenvalues \i of M. Moreover, 
Xm is a simple zero of the characteristic polynomial of M. 

It is known that a contact matrix is irreducible if and only if the neighbor graph 
it represents is strongly connected. A directed graph is called strongly connected if 
for any two vertices v i: Vj there exists a path starting at Vi and ending at Vj. Now 
we state and prove our main result of this section. 

Theorem 3.9. Let A G M 2 (Z) be an expanding similarity with characteristic poly- 
nomial f(x) = x 2 + px + q and T(A,V) be a disk-like CC tile. Then p(M) is the 
largest positive zero of the cubic polynomial 

% 3 ~ (\p\ ~ l)^ 2 — {\q\ ~ \p\) x ~ l<?l- 
Hence dim# (<9T) = 2 log p(M) / log |g| . 

Proof. Since | det(v4)| = \q\, it is more convenient to work with f(x) = x 2 ± px ± 
1 {p > 0, q > 2). Also we ignore those f{x) of the form f(x) = x 2 ± px — q (p > 
0, q > 2) as they cannot be characteristic polynomials of similarities (Lemma 13. 51) . 
In view of Appendix A or C, we can check that the contact matrix is irreducible if 
f(x) is one of x 2 ± 2x + 2, x 2 ± x + q, x 2 ± px + q (p > 2 excluding p = q = 2); 
and is reducible if f(x) = x 2 ± q. Hence we can apply Theorem 13.81 for all the cases 
except f(x) = x 2 ± q. 

Case (1) f(x) = x 2 + 2x + 2. The characteristic polynomial of the corresponding 
contact matrix M is (x — l)(x 2 + 2x + 2)(x 3 — x 2 — 2). We see that p(M) is equal 
to the positive zero of x 3 — x 2 — 2. By the cubic equation formula, we have 

, . [28 + (28 2 -l) 1 / 2 ] 1 / 3 + [28 -(28 2 -l) 1 / 2 ] 1 / 3 

p(M) = - — —± — — i — - — w 1.69562. 

3 
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Case (2) f(x) = x 2 — 2x + 2. The characteristic polynomial of the corresponding 
M is (x + l)(x 2 — 2x + 2)(x 3 — x 2 — 2). We have the same p(M) as in the case 
f(x) = x 2 + 2x + 2. 

Case (3) f(x) = x 2 + x + q. The characteristic polynomial of the corresponding 
M is (x - l)(x 2 + x + q)[x 3 - (q - l)x - q]. Notice that p(M) ^ 1. If p(M) = 1, 
then dim-^(<9T) = 0, which implies dT is totally disconnected (Proposition 2.5, |Fa] ) . 
This is not possible for the boundary of a topological disk. The zeros of x 2 + x + q 
are non-real. So p(M) is the largest positive real zero of a; 3 — (q — l)x — q. 

Case (4) f(x) = x 2 — x + q. The characteristic polynomial of the corresponding 
M is (x + l)(x 2 — x + g)[x 3 — (q — l)x — g]. The zeros of x 2 — x + q are non-real. So, 
as in the previous case, p(M) is the largest positive real zero of x 3 — (q — l)x — q. 

Case (5) f(x) = x 2 + px + q. The characteristic polynomial of the corresponding 
Mis (x— l)(x 2 +px + q)[x 3 — (p — l)x 2 — (q— p)x — q}. As explained before, p(M) ^ 1. 
The zeros of x 2 + px + q are either both negative or both non-real. Hence p(M) is 
the largest positive real zero of x 3 — (p — l)x 2 — (q — p)x — q. 

Case (6) f(x) = x 2 — px + q. The characteristic polynomial of the corresponding 
M is (x + l)(x 2 — px + g)[x 3 — (p — l)x 2 — (g — p)x — g]. We know that f(x) 
cannot have unequal real zeros. So we must have p 2 — 4g < 0. When p 2 — 4g < 0, 
the zeros of x 2 — px + g are non-real. Then p(M) is the largest positive real zero of 
x 3 — (p— l)x 2 — (q — p)x — q. When p 2 — 4g = 0, the two zeros of x 2 — px + q are equal. 
But the Perron-Frobenius eigenvalue should be a simple zero of the characteristic 
polynomial of M (Theorem 13 .8[) . so p(M) is also the largest positive real zero of 
x 3 — (p — l)x 2 — (g — p)x — q. 

Case (7) f(x) = x 2 + g. The contact matrix M is reducible. Its characteristic 
polynomial is (x 2 — q)(x 2 + q)(x — l)(x + l)(x 2 + 1). We see that p(M) = g 1 / 2 , which 
is also the largest positive zero of x 3 + x 2 — qx — q = (x 2 — q)(x + 1). 

Case (8) f(x) = x 2 — q. The contact matrix M is also reducible and its charac- 
teristic polynomial is found to be (x 2 — q) 2 {x + l)(x — l) 3 . As in the previous case, 
p(M) = g 1 / 2 , which is also the largest positive zero of x 3 + x 2 — qx — q. □ 

Remark 3.10. It is interesting to see that the signs of p and q do not matter in the 
calculation of dim h(8T) when A is a similarity. Notice also for the last two cases, 
f(x) = x 2 + q (\q\ > 2), we have p(M) = \q\ 1 ^ 2 . It follows that dim H {dT) = 1, as 



expected for the boundary of a parallelogram ( Proposition 1171) 



We observe that dim#(<9T) is independent of the choice of the vector v in the 
following sense. 

Corollary 3.11. Let A G M^i^) be an expanding similarity with characteristic 
polynomial f{x) = x 2 + px + q (|g| > 2). Let V = V(v, |g|) and V = V(v', |g|) be 
two CC digit sets such that each of {v, Av} and {v', Av'} is an independent set. If 
2\p\ < \q + 2\, then 

dim H (dT(A,V)) = dim H (dT(A,V)). 
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Proof. As 2\p\ < \q + 2\, both T(A,V) and T(A,V) are disk-like CC tiles (Theorem 
II. II) . Hence the corollary follows from Theorem 13.91 □ 



Remark 3.12. We conjecture that Theorem \3.9i and Corollary \ 3.11\ are also valid 
when 2\p\ > \q + 2\, i.e., T is non-disk-like. The major difficulty in justifying these 
conjectures is that, in general, there is no upper bound on the number of neighbors 
of a non-disk-like CC tile |DJNj . 



4. Appendix A: Neighbor Graphs 

Let f(x) = x 2 ± px ± q (p > 0, q > 2). The neighbor graphs of disk-like tiles are 
classified by f(x) and listed below. 




Figure 2. 




Figure 3. 
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(9-D 



(a) f(x) = x 2 + px + q, p > 2, 2p < q + 2 (b) f(x) = x 2 - px + q, p > 2, 2p < q + 2 
(excluding q = p = 2) (excluding q = p = 2) 

Figure 4. 



-G>+i), 



-p 



9-1 



" 

Av + pv 


q-p 




Av + ( 





fj+l »+■} 



-aft 













-(p+i) 



-p 



Av — pv 




-Av + pv 








1-V X 



q-p-1 



-p-l) 



Av - (fl. + l)w 


9 -P-l 


-A> + (p + l)u 




-w-p-i) 



(a) f(x) = x 2 + px - q, p > 1, 2p < q - 2 (b) f(x) = x 2 - px - q, p > 1, 2p < q - 2 

Figure 5. 




-i i 



/ 




-i i 



-V \ 1 



j4ll - 



— ife + 2d 



(a) f(x) = x 2 + 2x + 2 



(b) f(x) = x 2 -2.T + 2 



Figure 6. 
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5. Appendix B: Graph-directed Sets 

Let f(x) = x 2 ± px ± q (p > 0, q > 2). The graph-directed sets representing 
boundary dT are classified by f(x) and listed below. 

(1) f(x) = x 2 + q. Convention: u\ — v, U2 — Av — v, U3 — Av, M4 = Av + v. 



(2) f(x) = x 2 







9-1 q-1 

\J(T U2 +jv)U |J(T U3 
j=i 3=0 


9-2 

+ jv)u\J(T U4 +jv) 

3 = 


AT U2 




T 

1 -114 




AT U3 








AT U4 




T 




AT— U1 




9-2 9-1 

]J(T. U2 +jv)u]J(T. 

3=0 3=0 


9-1 

-n 3 +jv)u\J{T- U4 +jv) 
3 = 1 


AT- U2 




T Ui +(q- l)v 




AT— U3 




T U1 +(q- l)v 




AT— U4 




T-u 2 +(q- l)v 




— q. Convention: u\ — v, u<i - 


- Av — v, u 3 = Av, Ul 


AT U1 




9-1 9-1 

\J(T U2 +jv)u\J(T U3 

3=1 j=0 


9-2 

+ jv)u{J(T U4 +jv) 

3=0 


AT U2 




T- U2 +(q- l)v 




AT U3 




T U1 +(q- l)v 




AT U4 




T Ui +(q- l)v 




A r T— ul 




9-2 9-1 
\J(T-u 2 +jv)u\J(T- 

3=0 3=0 


9-1 

-u 3 +jv)U \J{T- U4 +jv) 

3 = 1 


AT- U2 




T 
± u 2 




AT- U3 




T— U1 




AT—U4 




T 





(3) f(x) = x 2 + x + q. Convention: u\ — v, u-i — Av, 113 = Av + v. 

9-1 9-2 

AT U1 = |J (T U2 + jv) U [J (T U3 + jv) 
3=0 j=o 

AT U3 T—ui 

9-1 9-1 

AT- U1 = [J (T_ U2 + jv) U |J (T_ U3 + jv) 

3=0 j=l 

AT- U2 = T U3 + {q- l)v 
AT- U3 = T U1 +(q- l)v 
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(4) f(x) = x 2 — x + q. Convention: U\ — v, u 2 — Av, u 3 = Av — v. 



AT U1 


9-1 q-2 

= \J(T U2 +jv)U \J(T U3 

3=0 3 = 


+ jv) 


AT U2 


T 

— -L U 3 




AT U3 


T 

— -t — ui 




AT— U1 


9-1 9-1 

= \J(T- U2 +jv)U \J(T. 

3=0 3 = 1 


-u 3 + jv) 


AT—u 2 


= T- U3 + (q- l)v 




AT— U3 


= T U1 +(q- l)v 





(5) f(x) = x 2 +px + q, p > 2, 2p < q + 2 (excluding p = q = 2). Convention:^ 
v, ii2 = Av + (p — l)v, u 3 = Av + pv. 

q-p 9-P-1 

AT U1 = \J(T U2 +jv)U |J (T U3 +jv) 

3=0 3=0 
p— 2 p—1 

AT U2 = (J (T- U2 + jv) U |J (T_„ 3 + jv) 

3=0 3=0 

AT U3 = T— U1 

9-1 9-1 

^_ U1 = U ( T -"2 +3") U U( T -"3+J^) 

3=P-1 3=P 
9-1 9-1 

AT- U2 = |J (T U2 +j«)U |J (T„,+jt;) 

3=9-p+l 3=1— P 

AT- U3 = T U1 + (g - l)w 



(6) /(x) = a; 2 — px + q, p > 2, 2p < q + 2 (excluding p = q = 2). Convention: 
Mi — v, U2 — Av — (p — l)v, u 3 = Av — pv. 

9-1 9-1 

AT U1 = (T U2 + jv)U \J(T U3 +jv) 

3=P-1 3=P 
P P—1 

AT U2 = |J (T U2 + j«) U |J (T„ 3 + jv) 

3=0 3 = 

AT U3 = T— U1 

q-p 9-P-1 

^T_ U1 = U(T_ U2 +jw)U |J (T_„,+jt;) 

3 = 3 = 

9-1 9-1 

^4T_ U2 = [J (T_ U2 +jw)U |J (T_„,+jt;) 

3=9-p+l 3=1— P 

AT- U3 = T U1 + (q - l)v 
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(7) f(x) = x 2 +px — q, p > 1, 2p < q — 2. Convention: U\ = v , u 2 = Av+pv, u 3 
Av + (p+ l)v. 

q-p-l q-p-2 

AT U1 = \J (T U2 +jv)U |J (T U3 +jv) 

3=0 j=0 

AT U2 = T U1 +(q- l)v 

q-l q-l 

AT U3 = |J (T U2 +jv)U |J (T„ 3 +ji;) 

3=9— P 3=1— P-l 

q-l q-l 

AT- Ul = (J( T -«2+^)u U ( T -" 3 +» 

j=p j=p+l 

AT— U 2 = Tu\ 

p-i p 
AT- U3 = |J (T- U2 + jv) U [J (T_ U3 + jw) 

3=0 3=0 



(8) /(x) = x 2 —px — q, p > 1, 2p < g — 2. Convention: «! = t>, w 2 = Av—pv, u 3 
Av — (p + l)v. 

q-l q-l 

AT U1 = \J(T U2 +jv)U |J (T„ 3 +jt;) 

j=p i=p+i 
^T U2 = T U1 + (g - l)i> 

q-l 9-1 

j=q—p 3=q—p-i 
q-p-l q-p-2 

3=0 j=0 

AT— ii 2 — T—ux 

p—l p— 2 

^_„ 3 = |J (T U2 + jv) U [J (T„ 3 + ju) 

3=0 3=0 



(9) /(x) = x 2 + 2a; + 2. Convention: U\ = v , u 2 = Av + t>, w 3 = At> + 2t>. 

AT U ^ — 

,4T U2 = T- U2 U T- U3 U (T_ U3 + v) 

AT U ^ — T— u -^ 
AT— Ul T—u 2 ~h f 

4T_ U2 = (T U2 +v)U T U3 U {T U3 + v) 
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(10) f{x) = x 2 — 2x + 2. Convention: U\ — v, u 2 — Av — w 3 = Av — 2v. 

AT U1 T U2 -+- v 

AT U2 = T U2 U T U3 U (T„ 3 + u) 



AT T 



^1T_ U2 = (T- U2 +v)U T- U3 U (T_ U3 + v) 



6. Appendix C: Contact Matrices 

Let f(x) = x 2 ± px ± q (p > 0, q > 2). The contact matrices (in table form) of 
disk-like tiles are classified by f(x) and listed below. 





V 


Av 


— V 


-Av 


Av — v 


—Av — v 


—Av + v 


Av + v 


V 





1 








3-1 








9-1 


Av 








1 

















—V 











q 





9-1 


9-1 





-Av 


1 























Av - v 

















1 








-Av - v 




















1 





—Av + v 























1 


Av + v 














1 












Table 4. f(x) = x 2 + q. 





V 


Av 


— V 


-Av 


Av — v 


—Av + v 


Av + v 


—Av — v 


V 





q 








9-1 





9-1 





Av 


1 























—V 











q 





9-1 





9-1 


-Av 








1 

















Av - v 

















1 








—Av + v 














1 











Av + v 




















1 





— Av — v 























1 



Table 5. f(x)=x 2 -q. 





V 


Av 


Av + v 


—v 


-Av 


—Av — v 


V 





9 


9-1 











Av 

















1 


Av + v 











1 








— V 














9 


9-1 


—Av 








1 











—Av — v 


1 


















Table 6. f(x) = x 2 + x + q. 
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V 


Av 


Av — v 


—V 


-Av 


—Av + v 


V 





q 


q-1 











Av 








1 











Av — v 











1 








—v 














q 


3-1 


—Av 

















1 


—Av + v 


1 


















Table 7. f(x) = x 2 - x + q. 





V 


Av + (p — l)v 


Av + pv 


—v 


—Av — (p — l)v 


— Av — pv 


V 





q-p+1 


q-p 











Av + (p — l)v 














p-1 


P 


Av + pv 











1 








—v 














q-p + 1 


q-p 


—Av — (p — l)v 





p-1 


P 











—Av — pv 


1 


















Table 8. f(x) = x 2 + px + q, p>2, 2p < q + 2 (excluding p = q = 2). 





u 


Av — (p — l)v 


Av — pv 


— V 


—Av + (p — l)u 


—Av + pv 


V 





q-p+1 


q-p 











Av — (p — l)v 





p-1 


P 











Aii — pu 











1 








— V 














q-p + 1 


q-p 


— /k> + (p — l)v 














p-1 


P 


—Av + pi; 


1 


















TABLE 9. f(x) = x 2 - px + q, p>2, 2p < q + 2 (excluding p — q — 2). 





V 


Av + pv 


Av + (p + l)v 


— 1! 


—Av — pv 


—Av — (p + l)f 


V 





q-p 


q-p-1 











Av + pv 


1 

















Av + (p + l)v 





P 


p+ 1 











—v 














q-p 


9 -P- 1 


—Av — pv 











1 








—Av — (p + l)v 














P 


p + 1 



Table 10. f(x) = x 2 + px - q, p>l, 2p < q - 2. 





V 


Av — pv 


Av — (p + l)v 


—v 


—Av + pv 


-Av + (p + l)u 


V 





q-p 


q-p-1 











Av — pv 


1 

















Av — (p + l)v 














p 


p + 1 


—v 














q-p 


q-p-i 


—Av + pv 











1 








-Av + (p + l)v 





P 


p+1 












Table 11. f(x) = x 2 - px - q, p > 1, 2p < q - 2. 
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V 


Av + v 


Av + 2v 


— V 


—Av — v 


—Av — 2v 


V 





1 














Av + v 














1 


2 


Av + 2d 











1 








— V 














1 





—Av — v 





1 


2 











-Av - 2v 


1 


















Table 12. f(x) =x 2 + 2x + 2. 





V 


Av — v 


Av - 2v 


— V 


—Av + v 


-Av + 2v 


V 





1 














Av — v 





1 


2 











Av - 2v 











1 








—v 














1 





—Av + v 














1 


2 


— Av + 2v 


1 


















Table 13. f(x) = x 2 -2x + 2. 
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